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ABSTRACT 

Spectral properties of super-Eddington accretion flows are investigated by means of 
a parallel line-of-sight calculation. The subjacent model, tak en from two-dimensional 
radiation hydrodynamic simulations bylOhsuga et alj l)2005f) . consists of a disc accre- 
tion region and an extended atmosphere with high velocity outflows. The non-gray 
radiative transfer equation is solved, including relativistic effects, by applying the flux 
limited diffusion approximation. 

The calculated spectrum is composed of a thermal, blackbody-like emission from 
the disc which depends sensitively on the inclination angle, and of high energy X-ray 
and gamma-ray emission from the atmosphere. We find mild beaming effects in the 
thermal radiation for small inclination angles. If we compare the face-on case with the 
edge-on case, the average photon energy is larger by a factor of ~ 1.7 due mainly to 
Doppler boosting, while the photon number density is larger by a factor of ~ 3.7 due 
mainly to anisotropic matter distribution around the central black hole. This gives an 
explanation for the observed X-ray temperatures of ULXs which are too high to be 
explained in the framework of intermediate-mass black holes. 

While the main features of the thermal spectral component are consistent with 
more detailed calculations of slim accretion discs, the atmosphere induces major 
changes in the high-energy part, which cannot be reproduced by existing models. 
We also conclude that, in order to interpret observational data properly, simple ap- 
proaches like the Eddington-Barbier approximation cannot be applied. 

Key words: super-Eddington accretion - spectral energy distribution - beaming 



1 INTRODUCTION 

With the constant improvement of observational techniques, 
more and more detailed information about accretion disc 
systems could be gained in the past decades. In particu- 
lar, a large number of ultraluminous X-ray sourc es (ULXs) 
have been discovered since the end of the 1980's iFabbiand 
[WS^)- These new astrophysical objects imposed a severe 
problem upon the existing general idea of accretion disc sys- 
tems: With a bolometric luminosity exceeding lO'^^ergs"^ 
(derived from X-ray observations), at least some of them 
show relatively low radiation temperatures (~ 0.1 keV). 
These systems have been suggested to be intermediate 
mass black hole (IMBH)^ sub-Eddington acc r etion disc sys- 
tems jMiller et alJl2003l: ICropper et aljr2004l:lRoberts et alJ 



l2005f) . However, from stellar formation theory, it is difficult 
to explain the formation of IMBHs. Therefore, many dif- 
ferent formation theories have been presented in the last 
years. A good overview can be found in lvan der Marell ll2004l) 



and iMiUer fc Colbertl i2004r) . Controversially, from time- 
variability observations, it seems to be likely that these ob- 
jects are instead low-mass X-ray binary systems \\An et alJ 
^2002)- Radio observations show that the distribution of 
ULXs can as well be fi tted by stellar mass black holes with 
mildly relativistic jets jKording et al.ll2004^ . 



With increasing observational data from X-ray satellites 
like ASCA, Chandra and XMM-Newton, large samples of 
ULX source s becam e avail ab le (e. g. IC olbcrt fc M us hotz kvl 
Jl999h: iM akishima ct~ aD (l200d) : iMiller et all ll20o3) : 
iKubotafc Makishima (ii200(t) ). They reveal that a distinct 
class of ULX sources exists, showing higher temperatures 
- sometimes e xceeding 1 keV - th a n can be ex plained by 
IMBHs (e.g. lOkada et al.l Jl998D : iMizuno et al. (i999J). 
Contrary, stellar mass Kerr black holes accreting above 
their Eddington lim it can account for these 
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JWatarai et al.l 1200 it lEbisawa et"ai] l2003l). Ahernatively , 
mild beaming could be important jKing et alj 1200 ll) . 

although the shape of the surrounding nebula indicates 
that it is illuminated by a nearly isotropic emission source 
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A controversial debate about the origin of 
ULX sources is still ongoing. 

Also in theory, a lot of investigations of accretion discs 
have been carried out in the past decades, not at last due 
to enhanced computational facilities. Systems hosting ro- 
tating black holes show higher temperatures and enhanced 
luminosities due to a smaller in ner boundary, compared to 
their non-rotating counter parts iGammie fc Pophamlll9S)3 
iPopham fc Gammielll99^ . To bypass the violation of the 
Eddington limit for stellar mass black hole ULXs, various 
mechanisms have been proposed. lAbramowicz ct al. (1988') 
exten ded the standard thin disc model ( Shakura fc Sunvaevi 
Il973l) by considering advective energy transport, leading to 
so called "slim discs". These discs allow for higher mass 
accretion rates, compared to thin discs with the same lu- 
minosity. The photon bubble instability mechanism permits 
higher accretion rates and high luminosities whi le preserving 
the global stability of the disc I^Begelmanll2002^ . An impor- 
tant extension of the slim disc mo del is the full treat ment 
of the photon trapping effects fsee lKat'3 (1977); Bcgcl manI 
l)l978h for the case of spherical accretion): For sufficiently 
high accretion rates, the gas inflow time becomes shorter 
than the photon diffusion time in the inner part of the disc; 
high-energy photons otherwise emerging from the inner disc 
region are partially swallowed by the black hole, leading to 
altered spectral energy distributi ons and lower luminosities , 
compared to a standard slim disc dOhsxiea et al.l2002Ll2nn3l . 
The question about the validity of the Eddington limit in 
accretion discs itself is widely discu ssed , yielding strong de- 
viations for these disc models (jlfeinzell ^r Duschl 2006'). 

Also the emerging disc spectra have been investigated 
widely in the past: Simple blackbody or modified black- 
body spect ra for supercritical accretion discs have been 
calculated jS zuszkicwicz ct al.' 'I99ef; 'Wa ng fc Zhol Il999l: 
lMineshige"et al.. ,200Q. : ..Watarai et al.. ,,200(ih . Furthermore, 
slim accretion disc spectra, including self-irradiation and 
self-occultation for self-similar solutions have been studied 
byfFukue ( 2000), while Watarai ct al. (2005) investigated 
the implications of geometrical effects an d general re l ativis- 
tic effects on the disc spectra. Further, iKawaguchil ll2003l) 
considered Comptonization effects in spectral calculations, 
finding significant spectral hardening occurring at large ac- 
cretion rates. In these approaches, however, anisotropy in ra- 
diation fields is not taken into account, although we naively 
expect mild beaming effects, i.e., radiation is likely to es- 
cape predominantly in the direction perpendicular to the 
disc plane. Moreover, the influence of the environment of the 
disc, e. g. its atmosphere, is not considered. It is well known 
that accretion discs own extended atmospheres, which in- 
deed have a strong influence on the emerging disc spectra. 
Photoionization of the accreti on disc surface by incid ent X- 
rays has been inve stigated bv iRevnolds et al.l (Il999h . while 
iDorrer et al.l lll996l) calculated disc spectra for a thin a-discs 
around a Kerr black hole, surrounded by a hydrogen atmo- 
sphere. 

In this study, we calculate the spectral energy dis- 
tribution of supercritical accretion flows based on the ra- 
diation h ydrodynamic (RHD) simulations computed by 
lOhsuga e t al. (200l). In Sect.|21 we describe the methods of 
calculations used in this investigation. Main aspects of the 
subjacent simulation data will also be briefly summarized 
there. We then present our spectral calculations in Sect.|21 




Figure 1. Sketch of the line-of-sight calculation 



Discussion will be given in Sect.|ll while Sect. |^ is devoted 
to conclusions and an outlook. 



2 MODEL SETUP 
2.1 RHD simulations 

In this study, we account for both a sophisticated disc model 
and the disc's atmosphere, computed in a self-consistent way 
within RHD simulations: we apply our calculations to the 
2D RHD simulation data from Ohsuga ct al. (2005). Start- 
ing with an empty disc and continuously injecting mass 
through the outer disc boundary, the authors simulated the 
structure of a supercritical accretion flow, until it reaches 
the quasi-steady state. The central object is given by a non- 
rotating stellar mass black h ole {M = lOMpj), genera ting a 
pseudo-Newtonian potential ijPaczvhski fc Wiitall98(l|^ . The 
viscosity is given by the classical a prescription. The mass 
input rate at the outer boundary (5007?g) strongly exceeds 
the Eddington limit, Mcxt = IOOOMb, where Me = Le/c^- 
The authors considered energy transport through radiation 
and advection and included relativistic effects in the radi- 
ation part. Note that photon trapping effects were auto- 
matically incorporated in the simulations. A gray computa- 
tion of the radiative transfer in the flux limited approxima- 
tion ijLevermore fc Pomranin3ll98ll ) was used. 

They found that the supercritical flow is composed of 
two parts: the disc region and the outflow regions above 
and below the disc. Within the disc region, the circular mo- 
tion as well as the patchy density structure are observed. 
The mass accretion rate decreases inwards (i. e. as matter 
accretes), roughly in proportion to the radius, and the re- 
maining part of the disc material leaves the disc to form an 
outflow. In particular, only 10% of the inflowing material fi- 
nally reach the inner boundary (3J?g), while the remaining 
90% get stuck in the dense, disc-like structure around the 
midplane or transform into moderately high-velocity out- 
fiows with wide opening angles. The outflows are accelerated 
up to ~ 0.1c via strong radiation pressure force. 

From the simulation data, the gas density p, its temper- 
ature Tgas and its velocity v are taken as input parameters, 
as well as the radiation energy density E. The methods of 
calculating other quantities, such as the radiation tempera- 
ture Trad, the source function Si, and the radiation pressure 
tensor P^, will be given in the following subsections. 
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2.2 Equation of radiative transfer 

Under the assumption of an observer being located at in- 
finite distance from tlie object, we calculate the emerging 
flux/luminosity as a function of the observers inclination 
angle O and azimuthal angle $. Here, O and $ refer to the 
spherical coordinate system which describes the computa- 
tional box/ 

More precisely, we adopt a parallel line-of-sight calcu- 
lation on a two-dimensional grid on the projected surface, 
seen by the observer (see Fig. for a better understand- 
ing). We start the line-of-sight calculation at a sufficiently 
high optical depth r^, start from the projected surface with 
initial intensity 1^=0 and with fixed direction cosine vector 

r= r(e,$). 

We focus on solving the radiative transfer equation nu- 
merically by a Runge-Kutta algorithm, considering the rel- 
ativistic corrections due to the high velocities of the gas. 
General relativistic effects such as gravitational lensing and 
gravitational redshift are not taken into account. We have 
to distinguish between two different coordinate systems: We 
describe the system in which the observer and the computa- 
tional area are in rest with Iv, I, etc., while the same quanti- 
ties are tagged with in the frame comoving with the local 
gas velocity v. Since the gas velocity strongly varies, the 
comoving frame depends on the position in the simulation 
box. 

In this framework, the relativistic equation of radiative 
transfer is given by 



+ 4'^;7^(^i5,„+Wo,P:f,j ^ (1) 
while, in the non-relativistic case, it reduces to 



where we assume the mass density of ions and electrons to 
be equal, p = pi = pe- For the scattering processes, we only 
consider electron scattering, given by 



sea ( P \ -1 

K = ot\ I cm 

Vmi, 



The relativistic transformation rules are 



h = — 

1^0 



i+{ c^^v-i-r] ^ 



(4) 

(5) 
(6) 
(7) 



with r being the Lorentz factor. 



2.3 Frequency-dependent radiation quantities 

Special attention is needed for deriving the quantities 
and El, (and therefore, by applying the flux limited dif- 
fusion approximation - see Sect. 12.41 - also P^). As the 
radiative transfer in the 2D RHD simulation is calculated 
in a gray approximation, the simulation data provides only 
frequency-integrated values for the radiation energy density. 
The matter distribution is described by the gas density and 
gas temperature. We assume local thermal equilibrium for 
the matter distribution and for the radiation field separately: 



(Tgas ) 



exp 



- 1 



Trad - ( f ) ^ 



a — radiation constant 



El, 



in 



Bua (Tr, 



(8) 

(9) 
(10) 



(/ ■ V)/. 



^abs Q 



(2) 



It is important to note that all quantities on the right hand 
side of Q are evaluated in the comoving frame, while those 
on the left hand side are given in the rest frame. In the 
equation of radiative transfer, Sv denotes the source func- 
tion for matter, Ei, the radiation energy density and P,y the 
radiation pressure tensor. Additionally, = i^t^^ + f^lf^- 

In this first approach, we restrict ourself to frequency- 
dependent absorption coefficients for free-free absorption 
processes and totally neglect bound-free absorption pro- 
cesses, = K^. This holds as a good approximation, 
since the gas temperature is mostly above 10^ K, and, hence, 
hydrogen is fully ionized. For simplicity, we do not con- 
sider for metal opa c ities. We adopt the formula given in 
iRvbick fc LightmanI ^7^, 



While © generally holds as a good approximation, Q and 
nul l have to be treated carefully: In a scattering dominated 
domain as it is the case in the underlying simulation data 
here, photons undergo multiple scattering and therefore ex- 
pand in space - accordingly, the radiation field is diluted and 
the photon number decreases. This implies that the average 
photon temperature will be underestimated by and IIUII . 
which will be discussed later (Sect. 01 . Therefore, we hence- 
forth focus our discussion on the relative changes of the fiux, 
photon energy and photon number due to the variations of 
the inclination angle O. 



2.4 Flux limited diffusion approximation 

To calculate the radiation pressure tensor, we apply the 
frequen cy-dependent flux limited diffu sion (FLD) approxi- 
mation (iLevermore fc Pomranin3ll981^ to the quantities in 
the comoving frame. In this context, the radiation pressure 
tensor can be expressed by 



kU' = 3.7 • io*r 



-1/2 /_P_ 



-hu/kT\ -1 

e ' cm 



(3) 



^ Note that the equatorial plane of the computational box is de- 
fined by the injection point of gas and its angular momentum 
vector. 



Pi/,0 — fi^ Ejjfl . 



(11) 



where fj, is called the Eddington-tensor. Its components are 
given by 

fj = i(l - U)S'^ + i(3/. - l)ntni . (12) 
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Here, denotes the normalized energy density gradient 



(13) 



Following iKlevI 11198911 . n1 and subsequent quantities can be 
expressed as functions of the energy density in the inertial 
frame. 

To close the resulting equations, the Eddington factor 
has to be determined. From the momentum equations, 
the relation between and \,j is given by 



fv — + X^TZ^, . 



\veA 



(14) 



The flux limitcr X^, itself can not be determined from the 
equations of radiative transfer, but has to be defined man- 
ually. In order to do so, two conditions have to be fulflUed. 
In the case of — > oo, the equations have to reduce to the 
classical diffusion limit, i. e. — » i. In the case of Xv ~* 0, 
the flux limiter must tend towards 1/TZv in order to ensure 
1^1,0 1 ^ cS,,o. 

Naturally, there exist multiple possibilities to describe 
the flux limiter X^. We adop t the common formulation 
from iLevermore fc Pomranina Hi 



A,(7^,) = 



2 + 7^^ 



37^^ + 7^2 



(15) 



2.5 Numerics 



In the present investigation we choose a single snapshot of 
the RHD simulation data after the simulation settled down 
in a quasi-steady structure. In this stadium, the structure 
does not change anymore in time in a signiflcant way, giving 
rise to consider our results as characteristic properties of 
such a system.^ 

The calculation of the disc spectra is performed as pre- 
sented above with and without relativistic corrections. As 
the simulation data is symmetric with respect to the az- 
imuthal angle we compute the spectra only in dependency 
of the inclination angle O. We investigate the results of the 
computation for different starting points (i. e., optical depths 
r^.start) for the line-of-sight calculation. For Ti,, start ~ 8, 
the results begin to saturate, leading to changes below 1% 
when starting at higher optical depths. We used r^, start = 10 
throughout all cases and validated the results with several 
integrations from Tj,, start = 10, yielding differences below 
10"^ 

The reason why such a low Ti,, start reveals the same re- 
sults than higher optical depths can be understood from the 
two extreme cases: 

(i) the gas is dense and cool with Tgas ~ Trad and small 
contributions of the gas to the total emissivity 

(u) it is diluted such that k'^'' <^ Kf^'' and the total emis- 
sion along the line of sight is completely dominated by the 
radiation fleld 

For all calculations it turns out that increasing gas temper- 
atures go hand in hand with dropping gas densities so that 
the opacity and total emission are governed by the radiation 
field. 

^ For details, we refer the reader to lOhsuga et all jSoO^. 



We divide the projected surface seen by the observer 
in a polar grid with coordinates (r, 93), see Fig. Both 
for the radial and the polar coordinate, we adopt a lin- 
ear grid with Nr = 100 and = 200 grid points. The 
discretization in frequency is taken to be logarithmic with 
Nv — 200 frequency values between 10^** Hz (0.5 eV) and 
10^^ Hz (50MeV). 



3 RESULTS 

3.1 Overall spectral properties 

Figure |5| shows the resulting spectrum uL^, for inclination 
angles O = 0, 7r/4 and n/2 with and without relativistic 
corrections. The luminosity is given by 



L.(e) = 47r / L{e,r,ifi)dA, 



(16) 



where A denotes the projected surface of the computational 
area, as it is seen by the observer. For low frequencies, the 
spectra only weakly depend on the viewing angle. Also, rel- 
ativistic corrections are unimportant for energies ^ 400 eV 
(u ^ 10^^ Hz). Contrary, for higher energies, the dependency 
on the viewing angle becomes stronger. For high inclinations, 
i. e. for an edge-on view on the system, relativistic correc- 
tions still remain unimportant, while they become drasti- 
cally visible for low inclinations, i. e. for a nearly face-on view 
on the disc. For both the relativistic and the non-relativistic 
case, an enhancement of the peak frequency and luminos- 
ity is observable for small inclination angles; although, this 
boost is much stronger when considering relativistic correc- 
tions. Furthermore, instead of a rapid drop of uL^ for high 
energies, a slower decline up to a plateau-like structure can 
be seen in all cases. 

Figure 13 illustrates again the spectrum for O = 0. To 
show that the observed plateau in Figs.|5|and|21is a result of 
the thermal emission ^'''^''5'^ from the hot gas in the photo- 
sphere, we calculate the spectrum without gas contribution. 



i. e. we set S,, 



■ everywhere. As it can be seen from the 



figure, the high-energy plateau disappears completely when 
neglecting the contribution from the hot gas, confirming our 
hypothesis. We want to remark that 

(i) this structure may by altered significantly if Compton 
scattering is taken into account, since this provides an effec- 
tive cooling mechanism for the gas. However, the inclusion 
of Compton scattering and the calculation of the decrease 
in Trad is beyond the scope of this work. 

(ii) it is unlikely to observe this plateau since the overall 
emission in this energy range is considerably low and the 
spectrum is dominated by the peak emission. 

To illustrate the influence of the disc's environment in the 
face-on case in Fig. 13 we calculate the spectrum for a 
"screened" central region (for R lOOiic, we set all phys- 
ical quantities to zero) and for the core region only (for 
R ^ 100-Rg, all physical quantities are set to zero). Here, 
R denotes the radial coordinate in the spherical coordinate 
system describing the computational box. The former case 
corresponds to a system where the inner 100i?G are entirely 
evacuated and so emission, absorption and scattering pro- 
cesses only exist outside the core region. Contrarily, the lat- 
ter case means that no emission, no absorption and no scat- 
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Figure 2. Disc spectra vLd for inclination angles = 0, 7r/4, 7r/2 with and without relativistic corrections 
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Figure 3. Disc spectra vL,^ for = with relativistic corrections. Beneath the normal spectrum (as in Fig. 0, we plot the spectra 
without gas contributions, for the core region only {R ^ IOORq) and for a screened inner region [R ^ lOOflc). Additionally, theoretical 
spectral shapes for thin and slim accretion discs are sketched 
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Figure 4. Total luminosity Ltot as a function of the viewing 
angle © for the relativistic and non-relativistic calculation 

tering takes place for R ^ 100-Rg, yielding the unaltered 
emission from the core region only. Obviously, the environ- 
ment of the disc has a rather strong influence on the emerg- 
ing spectrum. 

Furthermore, the theoretical spectral sh apes for a stan- 
dard thin a- disc dShakura fc Sunvaelll973t) and for a stan- 
dard slim disc jAbramowicz et al.l llQSai are indicated in 
Fig. 121 each time without consideration for self-irradiation, 
atmosphere, relativistic effects, etc. By just taking into ac- 
count the surface temperature TcS for these disc models and 
taking advantage of the face-on view (no self-occultation), 
vLi, oc u"^^^ for the thin disc case, and vL^ oc for the 
slim disc case (see, e. g.. iKato et alJ il998l Sect. 3.2.5)). 
These shapes do not coincide with our results, although they 
should be valid at least for the peak intensity region of the 
spectrum. 

3.2 Angular dependence of the luminosity 

In Fig. |1| we show the dependency of the total luminosity 

I'tot(e) = j L.{@)dv (17) 

on the viewing angle for both the relativistic and the non- 
relativistic calculation. The energy boost for small inclina- 
tion angles appears in both cases, although it is stronger for 
the relativistic calculation. Due to relativistic effects, the 
gain in luminosity compared to the non-relativistic calcula- 
tion LtotZ-t'toT""'' varies between 1.0 for 9 = 7r/2 and 1.9 for 
9 = 0. 

We thus conclude that Ltot is by a factor of ~ 6.4, at 
most, enhanced for a face-on observer, compared with an 
edge-on observer. The increase in total luminosity may be 
due either to an increase in photon number or an increase 
in average photon energy. Which one is more important? 

To answer this question and to outline the relativistic 
effects more explicitly, we plot in Fig. |^ the total photon 
number density ntot and the average photon energy IJu^) 
as a function of the inclination angle. From our SEDs, we 
calculate the photon number density by 

ntot(9)= j n,{Q)du = j ^^j^ dv (18) 

and from that, the average photon energy by (Jiv){Q) = 




{hu) 



0.2 U.4 0.6 O.S 1 1.2 1.4 



e 

Figure 5. Total photon number density ritot and average photon 
energy {hv) in the relativistic and non-relativistic case as func- 
tions of the inclination angle (see text for details) 



Table 1. Gain in total luminosity, photon number and average 
photon energy for the facc-on case (see text for details) 



e 




Ltot 




"tot 




{hv) 




rcl. 


non-rel. 


rel. 


non-rel. 


rel. 


non-rel. 





1.00 


1.00 


1.00 


1.00 


1.00 


1.00 


7r/6 


1.66 


1.26 


1.41 


1.17 


1.18 


1.08 


7r/4 


2.40 


1.59 


1.86 


1.39 


1.29 


1.14 


7r/3 


3.63 


2.14 


2.53 


1.72 


1.43 


1.24 


tt/2 


6.40 


3.40 


3.74 


2.35 


1.71 


1.45 



Ltot(9)/(cntot(9)). While relativistic effects become more 
or less unimportant in the the edge-on case, they cause an 
additional increase both in total number of photons orig- 
inating from the system and in average photon energy in 
the face-on case. Table summarizes the gain in total lu- 
minosity, photon number and average photon energy for the 
face-on case - this means, the ratio for each of the quan- 
tities for 9 = 0, compared to inclination angles between 
and 7r/2. 

The results given in TableQcan be explained physically 
as follows; Starting from the non-relativistic calculation, we 
find that 

• lower densities and therefore less effective absorption 
and scattering in the photosphere allow a deeper look into 
the hotter region for the face-on case, compared to the edge- 
on case. Hence, the average photon energy IJiv) is increased 
by a factor of 1.45. 

• photons can escape more easily through the diluted 
medium along the polar axis, while they get stuck in the 
dense disc-like structure concentrated in the midplane. The 
outflow is therefore collimated and the number of escaping 
photons raised by a factor of 2.35. 

At the same time, the (outflow) velocities of the gas close to 
the black hole {R ^ lOOiic) and around the polar axis are 
higher, which becomes important for the relativistic calcu- 
lation. 

• The frequency of the escaping photons is shifted from 
1/ to vq ^ by the relativistic Doppler effect, increas- 
ing the average photon energy additionally by a factor of 
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1.71/1.45 = 1.18, when comparing the face-on with the edge- 
on view. 

• Since the relativistic invariant is I/v^, the emerging in- 
tensity in the relativistic calculation, compared to the non- 
relativistic case is given by lo/I ~ (vo/u)^. One factor of 
vo/v directly goes into (hv) via the relativistic Doppler ef- 
fect, the remaining factor of {uq/i/)^ applies to the emerging 
photon number ~ Iv/{hv), raising it once more by factor 
of 3.74/2.35 = 1.59, when Q decreases from 7r/2 to 0. 

All in all, an observer located at O = 7r/2 only sees the 
emission from the outer part of the optically thick disc-like 
structure, which itself screens the relativistic effects in the 
inner region of the system. The radial velocities and also 
the azimuthal velocities are relatively low (v^ « 0.01c). For 
the mainly contributing part to the spectrum, the azimuthal 
velocity is (almost) perpendicular to the line of sight, there- 
fore the already weak relativistic effects are not detectable 
for an edge-on observer. In the face-on case, the highly rela- 
tivistic flow {va ^ 0.3c) can by observed due to the optically 
thin atmosphere above the disc. At the same time, the ra- 
dial velocity is pointing into the direction of the observer, 
leading to strong enhancements of the radiative flux at low 
inclinations. 



3.3 Blackbody fitting 

When spectral data of binary black hole sources are ob- 
tained, it is usual to flt them with blackbody (or disc black- 
body) spectra. We thus attempt a similar spectral fltting to 
our theoretically calculated spectra: We apply a non-linear 
least square flt to the emerging intensity by a black- 
body spectrum wit h temperature Tfit, altered b y a spectral 
hardening factor e llSoria fc PuchnarewicJ2002^ . The fltting 
function is then given by 



exp - 



(19) 



Note that the factor e * is introduced to ensure the same 
radiation energy loss: 



rp4 



(20) 



In order to account for stochastic fluctuations, we weigh the 
fltting coefficients by their relative intensity. So, the weight- 
function is given by 



w{v) = — , 

J tot 



Ii, dv . 



(21) 



Figure 1^ shows the results for the fitting temperatures Tat 
as a function of the inclination angle O. Additionally, we 
plot the surface-averaged radiation temperatures 



(22) 



at optical depths Ti, = 1 and Tv — 10. As mentioned earlier, 
the temperature of the radiation field will be underestimated 
by @ and llOII . We therefore concentrate on its relative 
changes for different inclinations and scale all temperatures 
to the fltting temperature Tfit at = n/2, where it is basi- 
cally the same for the relativistic and for the non-relativistic 
calculation. 




Figure 6. Fitted blackbody temperatures for the relativistic and 
non-relativistic calculation as a function of the viewing angle ©. 
Additionally, mean temperatures at = 1 and = 10 are 
shown. All temperatures are scaled by T-at at = 7r/2 



rel, Theia^O 

blackbody fit, rel, Theia=0 
rel, Thela=Pi/2 
blackbody fit, rel, Theia=Pi/2 
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Figure 7. Luminosity Li, and the corresponding blackbody fits 
for the face-on and edge-on case 



If neglecting relativistic corrections, the fltted black- 
body temperature is given roughly by the radiation tem- 
perature at a constant optical depth between 1 and 10. The 
blackbody temperature rises by a factor 1.3 when switch- 
ing from an edge-on to a face-on case. The spectrum is 
only weakly hardened compared to a Planck distribution 
at the same temperature Tfit'- The spectral hardening factor 
e adopts an almost constant value close to unity, e ~ 1.15, 
for all inclinations. 

When accounting for relativistic corrections, no surface 
of constant optical depth can be deflned any more: While the 
fltting temperatures resemble those of the non-relativistic 
case for high inclinations, they differ significantly for low 
inclinations, mirroring the above statement of stronger rela- 
tivistic effects for the face-on seen system. While the black- 
body temperature rises by a factor of 1.6, the hardening 
factor e stays almost constant around 1.15, like in the non- 
relativistic calculation. 

Finally, Fig. |7] shows the luminosity for the face- 
on and edge-on view and the corresponding blackbody fits 
in the relativistic case. Due to the weighting function, the 
peak intensity region is fitted quite well, while there are large 
deviations in the low-energy and high-energy part. 
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4 DISCUSSION 

The results presented in Sect.^lpermit to draw several con- 
clusions about the observational appearance of supercritical 
accretion disc systems. 

It is clearly not sufficient to consider only the disc and 
neglect its surroundings like its hot photosphere, outflow 
regions, etc. Their influence becomes most important in the 
high-energetic part of the spectrum {hv ^ 4keV). We find a 
plateau-like structure, independent of relativistic effects and 
of the viewing angle, which can be ascribed directly to the 
high gas temperature in the corona. Therefore, neither the 
basic thin disc spectrum, nor the basic slim disc spectrum 
fit to our results. 

Our results also confirm that the Eddington-Barbier ap- 
proximation, a common simplification of radiative transfer 
calculations for stellar atmospheres, cannot be applied in ac- 
cretion discs: In this approximation, one generally assumes 
that the emergent intensity along the line of sight is equal 
to the source function at constant optical depth r = 2/3. In 
our calculation, the main contribution to the emerging flux 
is produced at higher optical depths > 2/3; moreover, 
the exact value of t^s depends on the inclination angle. 

We observe an enhanced luminosity for more and more 
face-on seen systems, which is due to both enhanced aver- 
age photon energy and total photon number. Relativistic ef- 
fects alter the total photon number much more significantly 
(almost twice the non-relativistic treatment) than the aver- 
age photon energy. This can be identified as mild relativistic 
beaming. 

As mentioned above, due to ||5J and HUH , the result- 
ing fit temperatures underestimate the real temperature 
of the radiation field. Strictly speaking, the temperature 
of the radiation field should be determined at r* ~ 1 
(r* ~ \/ T-sca.j-abs'j g^jjj jjq^ at t « 1 (iu a scattering dom- 
inated domain, r ~ t^'^^). Hence, the radiation energy den- 
sity Ev will resemble more a shifted black body distribution 
with a colour-corrected temperature Tcoi, 



rather then HUH , where E^ 
ment of energy conservation. 



B.(rflt). 



(23) 

From the rcquirc- 



E 



■ du . 



(24) 



To get a rough idea on how much the derived temperatures 
are underestimated, we solve 1241 numerically for T^oi in the 
main emanating region of radiation (r ~ 10). This yields a 
correction factor ^ = Tcoi/Iflt ~ 10. With 



Tflt = [9.4 • lO'K ... 1.4 • lO^'K] (25) 

for O = [7r/2 . . . 0], this leads to colour-corrected tempera- 
tures in the range of 

Tcoi = [9.4 ■ lO'^K ... 1.4 • lO^K] (26) 

These temperatures would be consistent with the observed 
high temperatures of several ULX sources (Makish ima et alJ 
|200(J) that can not be explained in terms of intermediate 
mass black hole systems with sub-Eddington accretion rates. 
However, our approach is certainly too simplified to answer 
this "too hot accretion disc" puzzle in a satisfactory way. 

In this work, spectral hardening turns out to be neg- 
ligible. This may be due in parts to the assumption of 



le+()36 
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rel, non-fid, Theta^Pi/2 
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Figure 8. Disc spectra vLi, for © = 0, 7r/4, 7r/2, computed under 
the FLD approximation and complete isotropy 



Thomson scattering: Comptonizatio n effects are exp ected 
to harden the spectrum significantly zernv & Elvis 1987|; 
IRoss. Fabian fc Minestilge 1992: Kawaguchi..2003,'l . Then, if 
only the peak of the spectrum is observed, the absolute scale 
and therefore the spectral hardening factor e remains un- 
known and the observed temperature Tobs = eTcoi overes- 
timates the colour temperature Tcoi. Moreover, bulk mo- 
tion Compton scattering is known to alter photon ener- 
gies due to the angular red i stribution of the scattered pho - 
tons (iPsaltis fc Lambll997l) . ISocrates. Davis fc Blaej(l2004^ 
showed that turbulent Comptonization produces a signifi- 
cant contribution to the far-UV and X-ray emission of black 
hole accretion discs. 

A weak point in our investigation is the application 
of the flux limited diffusion approximation instead of solv- 
ing the full momentum equations: In this approximation, 
several terms in the equation of radiative transfer Q, like 
1/c^ {DF/Dt) with F being the absolute value of the flux, 
are dropped. These terms are of the order of v/c and may 
contribute to the relativistic effects we find in our spectral 
calculations. By calculating the emerging spectra under the 
classical diffusion limit (i. e. complete isotropy, A = 1/3), we 
find only little influence of the FLD approximation at all 
(see Fig. IHJ. Thus, the inconsistencies invoked by applying 
the flux limited diffusion approximation do not affect our 
results in a significant way. 



5 CONCLUSION AND OUTLOOK 

Our radiative transfer calculations, based on the 2D RHD 
simulation of highly accreting supercritical discs including 
the photon trapping mechanism, show that the interpreta- 
tion of observed disc spectra is not a straightforward task. 
Especially, we find moderate beaming effects when the sys- 
tem is viewed from nearly face-on; i.e., the average photon 
energy is larger by a factor of ~ 1.7 in the face-on case than 
in the edge-on case due mainly to the Doppler boosting. 
Likewise, the photon number density is larger by a factor of 
~ 3.7 because of the anisotropic matter distribution around 
the central black hole. Interpreting observations thus has to 
be done in a more sophisticated way than one may expect 
from basic disc models: It requires a careful treatment of the 
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radiative transfer with consideration for the discs surround- 
ings. 

We assume that both the gas and the radiation field 
separately stay in local thermal equilibrium. Although the 
weak coupling of matter and radiation (ftabs ^ Hscn) sup- 
ports this assumption, it remains questionable and also un- 
derestimates the temperature of the radiation fi eld. It is im- 
portant to note that previous invest i gation s by I Wang et alJ 
i|l99<i) : lFukii3 l|200(t : IWatarai et alJ l)200,'j ) also rely on this 
approximation; nevertheless, their results differ in a signif- 
icant way. Solving the crux of assuming LTE for the gas 
and for the matter distribution at the present stage is not 
possible, because it requires frequency-dependent RHD sim- 
ulations. 

In a next step, Compton scattering has to be included as 
well as frequency-dependent absorption for both continuum 
(bound-free absorption may become relevant in the low- 
energetic tail of the SED) and line processes: From the ob- 
servational side, emission lines, especially the K-shell tran- 
sitions of iron, are a prominent feature in accretion disc sys- 
tem and comprise m any details about the observed object 
(see iRevnolds et al.l iITqQQ ); RoYnolds. ( 200^ for example). 
Beneath the effects on the high-energetic part of the spec- 
trum mentioned before, Compton scattering will provide an 
efficient cooling mechan ism for the hot gas. 

As pointed out by IWatarai et alj ll2005h , also general 
relativistic effects should be considered in the vicinity of the 
black hole, which will primarily affect the spectra of face-on 
seen systems. 
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